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Abstract. The role of noise in the transport properties of quantum excitations is a 
topic of great importance in many fields, from organic semiconductors for technological 
applications to light-harvesting complexes in photosynthesis. In this paper we study 
a semi-classical model where a tight-binding Hamiltonian is fully coupled to an 
underlying spatially extended nonlinear chain of atoms. We show that the transport 
properties of a quantum excitation are subtly modulated by (i) the specific type (local 
vs non-local) of exciton-phonon coupling and by (ii) nonlinear effects of the underlying 
lattice. We report a non-monotonic dependence of the exciton diffusion coefficient 
on temperature, in agreement with earlier predictions, as a direct consequence of the 
lattice-induced fluctuations in the hopping rates due to long-wavelength vibrational 
modes. A standard measure of transport efficiency confirms that both nonlinearity 
in the underlying lattice and off-diagonal exciton-phonon coupling promote transport 
efficiency at high temperatures, preventing the Zeno-like quench observed in other 
models lacking an explicit noise-providing dynamical system. 
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1. Introduction 

Transport of quantum excitations in complex low-dimensional systems is a topic of 
paramount importance in many physical contexts, such as semiconductor nanowires 
and nanotubes [1] , metallic wires [2] , ultracold atom systems in one-dimensional optical 
lattices [3], quasi one-dimensional organic superconductors, such as semiconducting [4] 
and metal carbon nanotubes [5], 7r-conjugated polymers [6] and more complex quasi- 
1D nano-architectures for modern technological applications [7], including plastic light- 
emitting devices and organic solar cells [8]. 

Due to their crucial role in directing light energy to reaction centers during the 
early stages of photosynthesis [9], excitons occupy a prominent role among the studied 
quantum excitations [10]. The exciton concept was introduced in solid-state physics 
by Y. Frenkel in 1931 [11]. However, it was not until the 1948 seminal paper by A. S. 
Davydov [12] that this idea was applied to geometry-determined molecular systems, such 
as molecular crystals. These and related studies have paved the way for the investigation 
of exciton transport in light-harvesting bionrolecules, which contain embedded networks 
of light-absorbing pigments [10]. 

Considerable boost to the investigation of exciton transport in bionrolecules has 
been brought about by recent advances in 2D photon echo experiments, which have 
revealed unusually long decoherence times for excitons in light-harvesting complexes [13, 
14, 15,16] and conjugated polymer systems [17]. Moreover, theoretical evidence has been 
accumulated that noise in certain regimes could act as a protective factor for quantum 
coherences [18, 19, 20], increasing suitably defined measures of quantum efficiency 
related to the transfer of electronic excitation energy from a chromophore to a distant 
one. In turn, such findings have corroborated more detailed investigations of the coupled 
dynamics of exciton transfer and protein vibrations, pointing at a functional role of 
specific vibrational modes in promoting possibly function-related, long-lived electronic 
coherences [21, 22, 23, 24, 25, 26, 27, 28]. 

However, despite the great experimental and theoretical advances, many 
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fundamental questions remain open. In particular, pinpointing the structural 
determinants of the electron-phonon coupling that seem to provide noise-induced 
protection of exciton dephasing remains a challenging task. Moreover, few studies 
have addressed the role of the dynamical determinants of such mechanisms, i.e. 
the influence of specific inter-atomic and inter-molecular potential energy terms 
beyond the harmonic approximation. Importantly, nonlinearity is known to play 
a subtle role in many quantum transport processes, from heat conduction [29] and 
vibrational energy transfer [30, 31, 32] to photon-assisted electronic transport in different 
nanostructures [33]. More generally, it is well known that nonlinear effects modulate 
non-trivially transport in disordered dynamical systems [34, 35, 36, 37, 38]. 

Another issue of paramount importance that needs further investigation concerns 
the details of how the environment ( e.g. the degrees of freedom of the protein in 
photosynthetic complexes) couples to the quantum degrees of freedom. In a tight- 
binding perspective, where the quantum excitation is characterized by a given set 
of site energies {e*} and hopping rates {Ji\, this means investigating the effects 
of the specific functional dependence of {ej,j7i} on the degrees of freedom of the 
environment (dynamical disorder). For example, letting only site energies fluctuate with 
the environment amounts to adding a pure dephasing term to the unitary evolution 
in the Liouville equation for the time evolution of the one-particle density matrix of 
the quantum excitation. It has been proved that this leads to noise-enhanced transfer 
efficiency [19, 20]. However, it has been observed that pure dephasing is not a physically 
realistic scheme of coupling and that in general fluctuations of the hopping rates, even if 
smaller than those of the site energies, can have a considerable impact on the dynamics 
of quantum excitations [39, 40]. For example, this is the case of high-mobility organic 
crystals, such as pentacene and rnbrene, where large fluctuations in the hopping rates 
occur at room temperature [41, 42], In this case, it has been shown that Zeno effect at 
high dephasing rates is suppressed and one recovers asymptotic mobility of the quantum 
excitation at increasing noise, albeit with a diffusion coefficient that decreases with 
temperature [43, 44, 45]. 

Motivated by the above described open questions, in this paper we adopt a semi- 
classical modeling strategy to investigate the effect of noise on the mobility and transfer 
efficiency of a quantum excitation coupled to the vibrations of a one-dimensional atomic 
chain. More specifically, the two main points that we wish to address are: (i) the role 
of nonlinearity of the interatomic potentials of the underlying lattice and (ii) the role of 
local versus uon-local dynamical disorder. 

The structure of this paper is as follows. In section 2 we briefly discuss existing 
modeling strategies and present the microscopic model that we study in this paper, 
describing a quantum-mechanical quasiparticle ( e.g. an exciton or an electron) hopping 
on a one-dimensional lattice |. In section 3 we present the main results regarding 
the spreading properties of an initially localized state in a thermalized one-dimensional 

f Without loss of generality, we will refer to the quantum quasiparticle in the following as the exciton. 
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chain. In section 4 we analyze the quantum transport efficiency in our model in the 
presence of finite local recombination rates. Finally, in section 5 we summarize our 
conclusions and further discuss our main results. 

2. Exciton-lattice coupled dynamics 

The time evolution of non-isolated (open) quantum systems, i. e. systems in contact with 
a thermal bath, is an extremely difficult problem in general [46, 47]. Several approximate 
schemes have been proposed, including master equation approaches [48, 44, 49, 50] 
usually based on the projection operator technique [51], different non-perturbative 
methods [52, 53, 54] and path-integral based methods [55, 56, 57]. In some cases the 
exact solution of the master equation can be determined analytically [58]. 

Other molecular modeling approaches allow to consider a greater amount of 
microscopic detail through ab-initio simulations of both the quantum and the bath 
(e.g. the protein) degrees of freedom. These techniques combine molecular dynamics 
simulations for the dynamics of the environment with the time integration of the 
Schrodinger equation for the reduced quantum system, based on quantum electronic 
structure calculations [59, 60, 61, 62, 23]. Such methods, often referred to as 
Quantum Mechanics/Molecular Mechanics (QM/MM) are more sophisticated and 
detailed modeling schemes belonging to a more general family of modeling strategies, 
where the full quantum evolution of the system is parametrized by the classical 
coordinates of the underlying lattice/protein, that evolves in parallel according to 
Newton equations. Such schemes have been applied with success to a variety of 
problems, including energy and charge transport in polypeptide chains [63] and other 
biomolecules [31, 64, 65, 66]. 

In the tight-binding approximation, the most general Hamiltonian governing the 
propagation of an exciton coupled to a one-dimensional lattice can be written as 




) BlB m + ]T 


(1) 


m 


mn 


where M is the bead mass, B\ is an exciton creation operator at site n, and u n 
is the displacement of the n-th mass with respect to its equilibrium position. The 
term V (u n ,u m ) represents the total potential energy of the nodes m and n including 
possible onsite terms. The energies J nm (u n , u m ) describe the modulation of the quantum 
Hamiltonian due to the fluctuations of the underlying chain. In principle, depending on 
the physical nature of the system, both the site energies (n = m terms) and the hopping 
rates (n ^ m terms) can be be influenced by the vibrations of the chain atoms. 

The Hamiltonian in Eq. (1) describes a general class of semi-classical models, 
including the Davydov [67, 68], the Holstein [69, 70], and the Su-Schrieffer-Heeger 
(SSH) [71] Hamiltonians, that have been employed to describe the dynamics of different 
kinds of quantum excitations in a variety of physical systems. In such modeling schemes, 
one treats the exciton as a quantum mechanical particle, while describing the oscillations 
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of the lattice classically. Thus, the equations of motion (EOM) are given by a time- 
dependent effective Hamiltonian for the exciton, which depends on the lattice variables 
u n . The EOMs for the chain are those of a set of coupled oscillators driven by the 
exciton wave function. 

In this paper we specialize to a nearest-neighbour tight-binding scheme with fully 
fluctuating parameters (i.e. both site energies and hopping rates), coupled to a nonlinear 
Fermi-Pasta-Ulam (FPU) chain [72]. The FPU potential can be derived as a Taylor 
expansion of a generic nearest-neighbour interaction potential with respect to the 
equilibrium positions of the chain. However, for the sake of simplicity we neglect 
cubic terms, which are known to give rise to specific topological (kink) excitations in 
the system, where also the equilibrium position of atoms are shifted and also to more 
complex combined breather-kink modes [73]. 

The Hamiltonian governing the system is the sum of two terms H = H e + Hi. 
H e is the Hamiltonian for an exciton propagating on the lattice in a given dynamical 
configuration, while Hi is the lattice Hamiltonian. Explicitly, we have 

L L 

He = Ys en( u)BlB n + Y, Jniu) (i^+A + 5^n + l) (2) 

n=1 n=l 


and 


*■ = £ 


7^ + 2 ( U n+l ~ U n ) 2 + ^( u n+l 


n =0 



( 3 ) 


where L is the lattice size, M is the mass of the atoms and (3 is the anharnronicity 
parameter. The operator B n (T>f) destroys (creates) an exciton at the position occupied 
by the n-th bead. 

According to our scheme, the energies appearing in the exciton hamiltonian H e are 
renormalized over time by the lattice fluctuations, namely 


£n E n T Xe(^«+1 Un— l) i 


( 4 ) 


and 

Jn 'fn t Xj(^bi+1 'Un) ■ (5) 

In the above equations, the 0 superscript refers to the unperturbed values of site energies 
and hopping integrals, while the parameters Xe and Xj gauge the strength of the exciton- 
lattice coupling. In the following we will consider E® = J® = 1 V n, unless specifically 
stated otherwise. Notice that with the above choice of coupling the total Hamiltonian is 
invariant under global spatial translations of the u n , so that the total momentum of the 
lattice is conserved during the evolution. Moreover, the invariance of the Hamiltonian 
under phase transformations of the variables B n guarantees the conservation of the total 
excitonic probability. Eqs. (4) and (5) are a natural way to minimally couple the exciton 
to the chain, as they can be seen as the first order term of a Taylor expansion of the 
effective local energies and exchange integrals in powers of u n . 
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In our semi-classical approach, only the exciton state is treated quantum- 
mechanically, while the lattice variables u n evolve according to the laws of classical 
mechanics. As we work in the single-exciton manifold, in order to determine the time- 
evolution of the system we consider a trial wave-function for the exciton in the form 


bjj°) 

n 


( 6 ) 


and derive the appropriate coupled EOMs for the coefficients b n and lattice variables u n . 
The amplitudes b* n define the wavefunction in the basis of lattice sites. In the following 
we will always impose the normalization condition Y2 n \bn\' 2 = 1- The time evolution of 
the coefficients b n follows directly from Schrodinger’s equation, 

dibn 

^nbn Jn(b n -\-\ ~b 6 n _x) (7) 


where the parameters J n and e n are those of Eqs. (4) and (5), and depend on the lattice 
variables. The EOMs of the underlying lattice are obtained from the expectation value 
of the Hamiltonian on the exciton wavefunction, that is 


Mu n 


d 


du r 


{mw 


F + F e 

n 1 n 


( 8 ) 


with 

F n n(u n -|_i T u n — i 2 Un) T /3 [K +1 T {u n —i w n ) ] 

K = 2 X jRe(b* n+1 b n - b* n _-J) n ) + XE(b* n+1 b n+1 - &*_!&„_i) 

In the following we will take k — M = 1 so that the upper frequency of the linear 
spectrum of the chain is oj 0 = 2. Moreover, throughout this paper we consider periodic 
boundary conditions. It is worth noting that Eqs. (7) and (8) also correspond to the 
deterministic version of a set of equations that can be derived through a path-integral 
approach [55]. Furthermore, we observe that, despite our choice of a bilinear exciton- 
lattice coupling, the ensuing time evolution governed by Eqs. (7) and (8) is a nonlinear 
one. 

For the analyses reported in the following, we evolved numerically the coupled 
EOMs (7) and (8) starting from different initial conditions and referring both to 
equilibrium and non-equilibrium setups. For this we used a standard 4 — th order Runge- 
Kutta algorithm with a time step dt = 10 -3 . In particular, we prepared the lattice 
initial condition by sampling a characteristic configuration of the variables u n and p n 
representative of a finite temperature T. This task is accomplished by thermalizing the 
lattice via a Langevin heat bath [74] at temperature T for a sufficiently long transient 
time to- Explicitly, the Langevin thermalization is achieved by augmenting the free 
lattice EOMs with a suitable friction term and a stochastic force, namely 


Mu n = 


K ~ Wn + vErjnft) 


( 10 ) 
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where 7 defines the coupling strength of the reservoir and £ n (t) is a Gaussian white 
noise with zero mean and unit variance. Notice that the present setup corresponds to 
the thermalization of each site of the lattice chain with an independent heat bath at 
temperature T (measured setting the Boltzmann constant ks = 1)- Moreover, for all 
the simulations we will set t 0 = L and 7 = 1 . This choice of parameters guarantees 
an efficient thermalization of the lattice in the whole range of parameters studied in 
this paper. Finally, once the lattice has reached a stationary, thermalized state, the 
Langevin reservoir is disconnected in order to sample the microcanonical dynamics of 
the total (lattice + exciton) system. 

2.1. Equilibrium spectral analysis 

The explicit treatment of the lattice dynamics by means of Eq. ( 8 ) represents a simple 
and direct way to include more realistic, explicit spatio-temporal correlations in the 
environmental noise that perturbs the exciton evolution, as opposed to more abstract 
treatments where the environment only enters the picture as a spectral density. It 
is therefore interesting to provide a spectral characterization of the lattice dynamics, 
focusing on equilibrium stationary states at a given temperature T. The same analysis 
also allows us to explore the near-equilibrium exciton transfer processes. 

In order to show the relevant transport properties of the lattice, it is instructive to 
look at the power spectra of its long-wavelength Fourier modes. These are reported in 
the upper panel of Fig. 1. Simulations have been performed evolving the system in the 
presence of a Langevin thermal bath at temperature T [74] for a transient time to = L. 
We remark that, according to our prescriptions, time in our simulations is measured 
in units of \JM/n. The external heat bath was switched off at t = to and the power 
spectrum was then computed by sampling the microcanonical dynamics over an interval 
of 2 16 temporal units and averaging over different thermalized initial conditions. 

For a white-noise signal, one would obviously observe a flat spectrum. The sharp 
peaks that are visible in Fig. 1 flag the nontrivial propagation of correlations inside the 
nonlinear chain. Their presence is also closely related to the anomalous heat transport 
properties observed in the FPU chain [74] and clearly shows that the lattice dynamics 
cannot be approximated by a diffusive uncorrelated process like a pure-dephasing noise. 
We conclude that nonlinearity in the interatomic potentials of the underlying dynamical 
system couples the exciton to a noise possessing a complex structure. This is likely to 
be the case a fortiori for excitons propagating within complex fluctuating biomolecules. 
The necessity of including the proper correlations in the noise beyond pure-dephasing, 
possibly encoded in the underlying lattice structure, appears therefore important. 

We then move on to examine the lattice spectrum in the presence of non-vanishing 
exciton-lattice coupling, as this can inform on the back-action exerted by the exciton on 
the underlying dynamical system. In this situation, the thermalization process and the 
following free dynamics were performed by evolving the coupled equations (7) and ( 8 ) 
with an excitonic initial condition corresponding to a random-phase delocalized state. 
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Figure 1 . (Color online) Equilibrium power spectra S(k,u>) of the lattice dynamics 
(power spectrum of the time series '}2 n u n (t)e zkn ) for low-fc Fourier modes with 
k = 27 tto/L and m = {1,2,4,8} for a thermalized nonlinear chain with T = 1, /3 = 1 
and L = 256. The upper panel refers to an isolated lattice, while the lower panel shows 
the lattice spectrum in a regime of a strong coupling with the exciton. Each spectrum 
is computed by averaging over 50 independent realizations of the dynamics. 


The corresponding lattice spectra are shown in the lower panel of Fig. 1. We find 
that, even in the regime of strong coupling, the relevant features of the characteristic 
peaks are essentially unchanged. The only difference with the zero-coupling case is a 
slight deformation of the low-frequency and low-wavenumber region of the spectrum. 
We therefore conclude that our combined exciton-lattice model exhibits environmental 
correlations that generally survive for finite temperatures and coupling strengths. 

It is instructive to carry out a similar spectral analysis also for the excitonic degrees 
of freedom. In Fig. 2 we show the power spectrum S e (k, u>) of the exciton amplitude field 
\b n (t)\ 2 in the presence of a thermal background at finite temperature. The spectrum 
of each normal mode k is well fitted by a Lorentzian distribution, a manifest evidence 
that sufficiently strong perturbations resulting from the underlying lattice dynamics 
produce a diffusion-dominated exciton transport. It is important to observe that the 
spectral analysis reported here allows to explore the out-of-equilibrium properties of the 
exciton-lattice chain only perturbatively, i.e. in the spirit of linear-response theory. On 
the other hand, many realistic situations are characterized by strong out-of-equilibrium 
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Figure 2. (Color online) Power spectra S e (k,uj) of the exciton amplitudes ( i.e . power 
spectra of the time series J2 n \bn{t)\ 2 e lkn ) for an exciton-lattice coupling \e = = 1 

in a chain with T = 1, /3 = 1 and L = 256 (with k = 2irm/L). Each spectrum 
is computed by averaging over 50 independent realizations of the dynamics of the 
underlying lattice. The dashed black line is a power-law with exponent —2 and 
indicates the characteristic diffusive behavior. 


conditions. For example, this is the case of the propagation of photosynthetic excitons 
created in light-harvesting antenna complexes following the absorption of a photon. 
Such physical scenarios properly correspond to far-from-equilibrium initial condition 
that cannot be included in linear-response schemes. For this reason, in the next section 
we discuss in detail non-stationary exciton transport arising from spatially localized 
initial conditions. 

3. Exciton spreading on a chain at finite temperature 

In this section we study the spreading of an initially localized exciton wave-function 
interacting with the chain. The initial conditions for the lattice are taken by sampling 
random velocities from a Maxwell distribution at a given temperature T, which is one 
of the parameters in these simulations. The initial displacements are set to zero. The 
lattice is then evolved for a transient time t 0 = L in the presence of a Langevin heat 
bath at temperature T that interacts independently with each site of the chain [74], 
After this thermalization process, the external heat bath is disconnected and the lattice- 
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exciton interaction term in the Hamiltonian is switched on. Furthermore, we average 
the time evolution of the same initial condition for the exciton over many independent 
trajectories, each corresponding to different initial conditions of the chain sampled from 
the same thermal distribution. In our picture, the decoherence of the exciton wave- 
function is brought about by averaging over many independent realizations of the explicit 
noise (the lattice dynamics). 

For sufficiently low temperatures and couplings, one can argue that the exciton 
evolution should be only weakly perturbed by the environment. Indeed, the exciton is 
found to spread over the chain almost ballistically in this regime, with a slow loss 
of coherence. However, one may speculate that for larger values of couplings and 
temperatures the non-Markovian nature of the noise acting upon the exciton and 
the nonlinearity of the dynamics may play a fundamental role in modulating the 
spreading of an exciton. For instance, very large couplings typically result in the 
emergence of immobile self-trapped states of nonlinear origin. The net effect is that 
a substantial amount of the vibrational energy gets pinned around a handful of sites in 
the chain, producing a pinning potential where the exciton self-traps, making de facto 
impossible any kind of exciton transport. Self-trapped states are well-known in many 
problems studied with models belonging to the same class as ours, such as the Holstein 
polaron [70, 69] and discrete breathers (DB) in dilute Bose-Einstein condensates trapped 
in optical potentials [75, 76, 77, 78]. 

The situation is perhaps more interesting at intermediate couplings, where it is not 
clear a priori over what timescale the spreading is diffusive and what is the dependence 
of the exciton diffusion constant D on the lattice temperature. 

In order to gather information on the fraction of lattice sites that are significantly 
occupied during the time evolution of the system, we compute the participation ratio 
n, defined as 

EJftnOOl 2 1 1 

EnIMOl 4 EJfrnWI 4 

With this choice of normalization, it is easy to show that n = 0 for a completely 
localized exciton wavefunction, while one has n = L — 1 for a perfectly uniform state. 
Therefore, one should think at n as an effective length , measuring the spatial extent of 
the exciton wavefunction over the chain. As such, diffusive spreading would be flagged 
by a law of the type n(t) oc f 1 / 2 , while ballistic propagation would correspond to a linear 
dependence on time, n(f) oc t. 

In Fig. 3 we show the typical time evolution of the participation ratio for different 
coupling strengths and T = 0.1. The first stage is a short transient {t < 1), where n 
grows quadratically in time. This super-ballistic evolution is characteristic of the very 
first stage of the time evolution of an initially localized wavefunction. This can be easily 
proved by writing down Eqs. (7) for an exciton initially sitting entirely at site n in 
a chain with J n = J V n. In this case it is not difficult to show that, if one defines 
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Figure 3. (Color online) Evolution of II(i) for T = 0.1 and different couplings in a 
chain with L = 1000 lattice sites. The dashed line is a plot of the initial super-ballistic 
propagation II(t) = (t/r) 2 with r = h/(2J). Dotted and dot-dashed lines indicate a 
power-law growth of the kind t a with a = 1 (ballistic) and 0.5 (diffusive), respectively. 
The inset shows the final profile of \b n \ 2 for the largest coupling \e = XJ = 10- A 
self-trapped state is clearly observable around the center of the chain. 


Q(t) = \b n (t)\ 2 - \b n ±i(t)\ 2 (with <2(0) = 1), then 


Q(t) = exp 



sin 2A cp A </> 1 dt' 


( 12 ) 


where r = h/(2J) and A0 = (f> n ±i — (f) n is the phase difference between the initially 
excited site and its neighbors. Without loss of generality, in the spirit of a Taylor 
expansion, we can assume A <f> ~ t/r in the early stages of the propagation. It is then 
straightforward to show that this immediately leads to a super-ballistic trend, namely 
11(f) ~ (t/r) 2 . It is apparent from Fig. 3 that this prediction is in excellent agreement 
with the simulations for low and moderate couplings. For larger values of the coupling, 
the physical law does not appear to change, while the time constant turns out to be 
renormalized. For example, for Xe = XJ = 10 we get 11(f) ~ (t/r') 2 with t' = r/3, which 
means faster super-ballistic propagation for large coupling strengths. It is intriguing to 
observe that this first super-ballistic stage could be physically relevant in many contexts. 
For example, in light-harvesting complexes one typically has J ~ 100 cm -1 , i.e. r ~ 170 
fs, which is of the same order of magnitude as the observed lifetime of quantum beats 
in 2 D photon echo spectroscopy at room temperature [13]. 

In the subsequent evolution for times greater than r = h/(2J), we can single out 
two main regimes corresponding to different dynamical situations. At low couplings, we 
recover the expected almost unperturbed evolution associated with ballistic spreading 
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of the exciton, i.e. II (t) oc t. The plateau observed at long times in Fig. 3 for 
Xe = Xj — 0-01 and 0.1 simply signals that the exciton wavefunction has reached 
complete delocalization (the lattice is finite) and no further spreading is thus possible. 
For intermediate couplings, it is clearly possible to observe a transition from a ballistic 
to a diffusive (II(t) oc f 1 / 2 ) regime. Overall, we conclude that for intermediate couplings 
one should always expect a ballistic-to-diffusive crossover, the time scale associated with 
it decreasing with increasing coupling. 

The situation changes rather dramatically for large values of the coupling strengths. 
In this case, the system enters a strongly nonlinear regime, where the localized initial 
condition triggers the spontaneous creation of a stable self-trapped state of nonlinear 
origin, which results in a very slow sub-diffusive transport (see also the inset of Fig. 3). 
Increasing the coupling further causes the initial amplitude of the exciton to stay 
permanently stored in a localized, time-periodic excitation of the system which is 
virtually decoupled from the rest of the system (cf. Ref. [77] for a comprehensive 
review on discrete breathers in nonlinear lattices). It is important to point out that 
the asymptotic stability of the self-trapped state depends both on the strength of the 
coupling xe — Xj and on the temperature T. In particular, for certain critical values 
Xc an d T c , such localized structures become unstable and get quickly destroyed by the 
thermal fluctuations of the lattice [79, 80]. On the other hand, the ballistic and diffusive 
regimes shown in Fig. 3 are not separated by a true dynamical transition. In fact, one 
can argue that for sufficiently long chains and times, any arbitrarily small interaction 
with the lattice will eventually cause a diffusion of the excitonic wavepacket. 

For the exciton nonequilibrium evolution reported in Fig. 3 we have also monitored 
the lattice kinetic temperature T\ measured after the Langevin heat bath has been 
disconnected. T k is defined as 

nm = 1[<pS(i)> - (pjrf] (i3) 

n= 1 

where the symbol (•) refers to a classical average over the set of independent lattice 
trajectories. Fig. 4 illustrates the evolution of T k for the same setup of Fig. 3, which 
corresponds to a temperature of the Langevin bath T = 0.1. Interestingly, in a wide 
region of coupling values that keep the overall system out of the strongly nonlinear 
regime, T k remains close to the Langevin temperature during the whole exciton evolu¬ 
tion. Conversely, the emergence of a stable discrete breather for large coupling strengths, 
Xe = Xj — 10 ( c f- blue diamonds of Fig. 3), produces a clear increase of T kj which is 
associated with the conversion of a substantial portion of lattice energy into (negative) 
exciton-lattice interaction energy. Accordingly, by virtue of the conservation of the total 
energy of the system, such energy transfer causes the lattice to heat up. However, this 
should be regarded as a finite-size effect. In general, we expect that, upon increasing 
the lattice length L, the heating effect becomes less and less important until it should 
eventually disappear in the thermodynamic limit L —>■ oo, since the breather interaction 
energy is localized over a finite number of lattice sites, whereas the lattice energy scales 
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Figure 4. Evolution of the lattice temperature T j. (see Eq. 13) during the exciton 
spreading dynamics shown in Fig. 3. 


linearly with L. Altogether, the above analysis confirms the consistency of the lattice 
dynamics as a well defined explicit thermal environment. 


3.1. The effective diffusion constant 

From the above discussion it is clear that in the intermediate coupling regime the 
asymptotic dynamics is diffusive. We now turn to analyzing in detail the properties 
of the diffusive spreading by a characterization of the diffusion constant D , defined as 

D = lim r 1/2 n(t) (14) 

£—>•00 

In Fig. 5 we compare the growth of fl(f) for increasing temperatures T and fixed values 
of the couplings Xe = Xj — 1 bi the intermediate regime. Interestingly, we find a 
nonmonotonic behavior for the dependence of the diffusion constant on the temperature 
T. More precisely, we observe a minimum located around T = 10, indicating the 
presence of a slowed-down spreading dynamics at intermediate temperatures. 

This phenomenon can be illustrated more effectively by measuring the average 
time t s it takes for the participation ratio to reach a certain threshold value fl s . In 
this temporal representation a minimum in the diffusion constant D corresponds to a 
maximum of the time t s . The value of fl s needs to be chosen in such a way as to 
avoid both the transient dynamics (typical of short times) and the saturation of fl(f) 
due to the finite size of the considered lattices (see again Fig. 5). Accordingly, we have 
chosen fl s = 100 for a chain of L = 1000 sites. The dependence of t s on temperature 
is illustrated in Fig. 6 for increasing couplings. For the lowest coupling considered, t s 
displays an initial growing stage, followed by a decrease at high temperatures after a 
maximum, which becomes more and more distinguishable upon increasing the coupling 
(note the logarithmic scale on the y axis in Fig. 6). Interestingly, this maximum appears 
to move towards smaller and smaller temperatures at increasing coupling strengths. This 
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Figure 5. (Color online) Time evolution of II(t) for \e = Xj — 1 and different 
temperatures of the lattice. The inset shows the behavior of the diffusion constant 
measured by fitting a power law to the asymptotic portion of the spreading. The 
dashed line is a plot of a power law with exponent 1/2. 


feature should be compared with the data displayed in the inset in Fig. 5, illustrating 
a minimum of the transport coefficient associated with exciton transport. A stationary 
point at the same value of temperature T & 10 is indeed recovered in both mobility 
indicators, D and t s . 

The sudden growth of t s in the low-temperature region for x — 2.5 flags the 
presence of self-trapped, breather-like excitations, which pin energy locally and thus 
slow down the relaxation process. This kind of localized states, however, are not present 
at the temperatures characterizing the maximum of t s , since the strength of thermal 
fluctuations is too large to sustain coherent localized nonlinear vibrations [79, 80]. 

The above analysis suggests that the nonmonotonicity displayed by transport 
coefficients with temperature is not related directly to nonlinear localization phenomena 
that pin energy down at high temperatures. It is therefore highly likely that the observed 
transport behavior is directly linked to the specific way the lattice vibrations couple to 
the parameters entering the exciton Hamiltonian. In fact, in the limit case of spatially 
uncorrelated white noise, it has been long known that local and nonlocal perturbations of 
the coherent Schrodinger dynamics can produce dramatically different behaviors for the 
spreading of an initially localized excitation [43]. In particular, it is known that purely 
local noise (in our picture, dynamical modulation of the site energies only) results in 
suppression of transport for large dephasing rates (a phenomenon often considered as 
an instance of the quantum Zeno effect [81, 20]). 

With these ideas in mind, we turned to examine the role of the effective noise acting 
on the free exciton dynamics as a consequence of the lattice thermal fluctuations. In the 
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Figure 6. (color online) The average time t s needed to attain a participation ratio 
n s = 100 as a function of temperature T for different values of the coupling strengths. 
The chain length is L = 1000. 


same spirit as the analysis performed in Ref. [43], we simplified the coupling between 
the exciton and the lattice by studying separately diagonal (involving site energies) and 
off-diagonal (involving hopping rates) interactions. 

In Fig. 7 we compare the spreading diffusion constant of a system exhibiting only 
diagonal coupling (xj = 0 ) with the one corresponding to a pure off-diagonal coupling 
(Xe = 0). Interestingly, the nonmonotonic behavior of D is present only in the latter 
case, while in the former we observe a monotonic decrease of D with temperature. This 
is precisely what happens in a master equation description a la Haken and Strobl [44] 
with nearest-neighbor Coulomb coupling, where 


D 


2yi a 2 + 


a 2 J 2 

ft 2 (7o + 37 i) 


(15) 


Here y 0 and 71 are the diagonal (pure dephasing) and off-diagonal noise strengths, J is 
the Coulomb hopping integral and a is the lattice spacing. The results of our simulations 
performed with xj — 0 are i n agreement with the prediction (15) with Xi = 0, that is, 
a value of D which decreases monotonically with temperature (70 in the language of 
Refs. [44] and [20]. In fact, the effect of purely local perturbations of the excitonic energy 
is to produce a diffusion of the Schrodinger phases that inhibits quantum transport. 
Eventually, in the limit of infinite interaction (i.e. infinite temperature), the quantum 
system remains frozen in the initial condition (D = 0) as a consequence of the complete 
randomization of the phases. 

The situation changes when we let the hopping rates be modulated by the lattice 
dynamics (xj 7 ^ 0, i.e. 71 7 ^ 0). As can be appreciated from Fig. 6 , the diffusion constant 
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Figure 7. (Color online) Temperature dependence of the spreading diffusion constants 
D in the presence of purely diagonal interaction (black circles) and purely off-diagonal 
interaction (red squares). Blue diamonds refer to the full model, with \e = 1 and 
Xj = 1. Other parameters are L = 1000 and /? = 1. 


displays a minimum, in agreement with the general prediction (15) and the transport 
becomes faster at increasing temperatures. In fact, since both phase and amplitude 
perturbations are now allowed for, the limit of infinite temperature corresponds to an 
arbitrarily large diffusion constant. This result can be interpreted as a recovery of 
classical amplitude diffusion in the infinite temperature limit. On more formal grounds, 
it can be shown that such classical reduction allows one to reduce a generalized Lindblad 
equation [43] for the exciton density matrix to a classical Fokker-Planck equation for the 
amplitude probability distribution. Although the non-Markovian nature of the effective 
noise and the explicit nonlinearity make the analytical calculation of D(T,xe,Xj) 
extremely difficult, the numerical results reported in Fig. 7 clearly indicate that classical 
diffusion is recovered also when one considers the full model with both local and nonlocal 
couplings. 

Interestingly, we remark that at low temperatures pure-diagonal noise allows for the 
fastest transport. However, as non-diagonal noise results in a minimum, the situation 
reverses beyond a characteristic temperature when \j is switched on. In this case, 
the model with fluctuating coupling strengths becomes the one affording more rapid 
spreading at high temperatures. 

It is interesting to note that extended vibrational modes are required in order 
to observe a nonmonotonic behavior of the diffusion constant, signaling a nontrivial 
coupling between exciton spreading and collective modes of the underlying lattice. This 
can be appreciated by comparing our analysis with the results presented in Ref. [45], 
where the hopping rates in the tight-binding exciton Hamiltonian are modulated by the 
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Figure 8. (Color online) Temperature dependence of the spreading diffusion coefficient 
D for an harmonic chain (/? = 0) with pure diagonal coupling (xe = 1 and \J = 0). 
The red line is a power law fit with a function D oc T -7 with 7 = 0.58. 


dynamics of a set of uncoupled harmonic oscillators. In this case, the authors report 
values of the diffusion coefficient that decrease monotonically with temperature. This 
is possibly a consequence of the absence of coupling between the oscillators providing 
the noise. Alternatively, they might be just exploring the low-temperature regime, as 
defined by their parameters. 

From Fig. 7 one can also argue that the effective interaction experienced by the 
exciton for finite temperatures can not be mapped onto a spatially and temporally 
uncorrelated noise as in [43]. Specifically, we find that the spreading problem for pure 
local exciton-lattice interactions is ruled by a nontrivial power-law decay D(T) ~ T -7 
with 7 fa 0.3 (data in Fig. 7), whereas 7 = 1 for white noise [43]. In the absence 
of nonlinearity in the lattice Hamiltonian (/3 = 0), the characteristic exponent is 
found to be close to 7 ~ 0.6, as shown in Fig. 8. We therefore identify two different 
sources of slowing down in the transport, that are characteristic of the explicit lattice 
dynamics. The first one is due to the presence of spatio-temporal correlations in the 
lattice system. The second one is associated with explicit nonlinear terms in the lattice 
pairwise potential energies V (u n ,u m ). 

here put some additional consideration on the effect of the /3 nonlinearity. 

Overall, the above detailed analysis allows us to conclude that the high-temperature 
limit of excitonic systems interacting with noisy environments is crucially determined 
by the specific properties of the exciton-lattice coupling and by the nature of the 
spatio-temporal correlations that characterize the noise-providing underlying lattice. 
In particular, the presence of non-diagonal coupling is sufficient to suppress localization 
at high temperatures, which can be regarded as the semiclassical counterpart of the 
quantum Zeno effect observed in quantum master equation approaches [?]. 

An interesting consequence of the above reasoning is that the standard scenario for 
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noise-assisted quantum transfer efficiency [19, 20] may display novel structures when 
passing from local pure dephasing noise to more realistic models of coupling including 
amplitude-affecting terms in the non-Hermitian part of the Hamiltonian. We thus 
turn now to discussing the implications of our explicit-noise approach for the efficiency 
problem. 

4. Exciton energy transport efficiency 

A quantum excitation such an exciton has an intrinsic lifetime, which is dictated by the 
recombination rate q r associated with the specificities of the environment. For example, 
in light-harvesting systems y r is estimated to be about 1 exciton per nanosecond [20]. 
The quantum excitation is therefore damped as it spreads through the system following 
its excitation. It is interesting to provide a measure of efficiency associated with the 
transport of an exciton. This can be done by requiring that a sink exists at some specific 
location in the system ( e.g . allowing the excitation to be transferred to a neighboring 
equivalent system) and evaluating the probability that the exciton exits through the 
sink rather than decaying non-specifically due to recombination mechanisms. 

In a master equation description, a recombination probability and a sink appear 
as non-Hermitian terms in the time-evolution operator for the exciton density matrix. 
Similarly, in our approach we ought to add damping terms to the equations of motion (7). 
More precisely, we consider a chain where the site k is identified as a sink, characterized 
by a trapping rate T y r . Therefore, the modified equations of motion read 


ih Jnipn+1 bn—l) Hlff r T 5nfcT)6 n 


(16) 


Along the lines of previous studies, such as Refs. [20], [82] and [83], we use eqs. (16) to 
investigate the competition between the two mechanisms of exciton destruction, namely 
generic recombination (y r ) and exit through a specific channel (T). The whole idea 
is that an efficient transport is maximally effective in channeling the exciton rapidly 
through the specified exit site against the generic degradation due to recombination. 
For example, this might reflect an exciton leaving a light-harvesting complex through 
a specific pigment connected to the reaction center. Moreover, as done in the above- 
cited studies, we take the unperturbed exciton site energies in Eq. (16) as random, 
which has the well-known effect of inducing spatial localization of the exciton at zero 
temperature. 

The incorporation of non-Hermitian terms implies a loss of norm as time passes, 
until eventually the norm of the exciton wave-function reaches zero for infinite times. 
Accordingly, a measure of transport efficiency can be computed in the following fashion 



(17) 


One can easily prove that ij is finite and takes values in the interval [0; 1] as a 
consequence of the conservation of the total excitonic amplitude when = T = 0. The 
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Figure 9. (Color online) Exciton transport efficiency versus temperature in different 
kinds of lattices of length L = 200. The recombination constant and the sink 
dissipation are, respectively, 7 r = 10~ 4 and T = 10 _1 . The local unperturbed exciton 
energies E ° are randomly chosen from a uniform distribution in the interval [—0.5 ; 0.5]. 
The efficiency ry is obtained by averaging over 20 independent realizations of the static 
disorder on the site energies E® and of the dynamical lattice noise. 


integral in Eq. 17 can be computed numerically to any desired accuracy, as controlled 
by the amount of total exciton norm N(t) = ^7 |fcj(t)| 2 left in the system at time t. In 
our simulations we integrated equations (16) until the survival probability reached the 
threshold value of 10 -5 . 

In Fig. 9 we compare the transfer efficiency of our chain model with the 
one corresponding to a quantum dynamics in the presence of pure-dephasing (PD) 
noise [43, 20]. The lattice is initially thermalized at temperature T with the same 
procedure described in the introduction of Section 3, while the PD dynamics describes 
the interaction of the exciton with a classical incoherent external field (no explicit lattice 
in this case) that induces decoherence on the quantum system. This is accomplished by 
specializing Eqs. (4) and (5) to 

e n = E® + \/2T£ n (t), J n = (18) 

where C, n (t) is a Gaussian white noise with zero mean and unit variance and T effectively 
accounts for the dephasing rate of the process. In both cases the exciton is initially 
injected at one side of the chain, while a trap is located at the opposite end. Interestingly, 
in the low-temperature region the two systems display qualitatively the same behavior, 
namely a disruption of disorder-induced localization due to increasing dephasing (PD 
model), i.e. increasing thermal fluctuations of the underlying lattice in our scheme. 
However, at higher temperatures the efficiency of transport in the presence of explicit 
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(i.e. produced by the lattice) noise turns out to depend on the details of the underlying 
lattice dynamics. The first remarkable result is that nonlinear correlations in the lattice 
dynamics boost the efficiency in the Zeno-effect region (see diamonds versus squares 
in Fig. 9). This finding, although for somewhat different physical reasons, is in line 
with recent results where the importance of spatial correlations was demonstrated for 
transport in a PD-like model [83]. Furthermore, we note that such a behavior is 
consistent with the different scaling exponents of the diffusion constant D discussed 
in Section 3. 

The second important finding is that non-diagonal noise (i.e. non-zero coupling 
between the hopping rates in the exciton Hamiltonian and the lattice) suppresses the 
Zeno drop in the efficiency at high temperatures. This is in good agreement with the 
prediction of eq. (15), namely that the diffusion coefficient should be non-monotonic 
with temperature (i.e. non-diagonal noise, 71 ). We recall that in this regime, the 
master equation for the density matrix turns into a classical diffusion equation, which 
should guarantee diffusive transport [48], albeit possibly with a diffusion coefficient that 
decreases with temperature [45]. 

5. Conclusions 

In this paper we have studied a model describing the dynamics of a quantum excitation 
that propagates in a system at finite temperature. In our scheme, the quantum 
evolution is dictated by a tight-binding (TB) Hamiltonian, whose matrix elements 
are functions of the classical coordinates of an underlying one-dimensional lattice. 
We refer to this setting as an open quantum system with an explicit environment, 
which provides a direct source of noise to the quantum excitation, endowed with 
specific spatio-temporal correlations. While the main ideas behind this modeling 
scheme are not new (see for example Ref. [67]), our study contains important elements 
of novelty. Notably, we explicitly focussed on spotlighting the signatures of non¬ 
trivial spatio-temporal correlations of nonlinear origin expressed by the underlying 
lattice. Furthermore, we conducted an original investigation of the relative role of 
diagonal and non-diagonal exciton-phonon couplings. In both cases, we uncovered a 
rich phenomenology, prompting new directions of investigation. 

We first examined the spreading of an initially localized excitation. Our results 
show that the very first stage of the propagation is faster than ballistic, up to a time 
of the order of h/(2J), J being the magnitude of the hopping integrals in the TB 
Hamiltonian. The subsequent time evolution is characterized by a transition to a ballistic 
stage followed by a crossover to an asymptotic diffusive regime, which appears at earlier 
and earlier times as the exciton-phonon coupling strength is increased. However, for 
large values of the coupling the picture changes dramatically, as a self-trapped state of 
nonlinear character sets in after the first super-ballistic spreading. The result is that 
the transport is completely suppressed in this regime as a sheer nonlinear effect. 

An analysis of the diffusion coefficient D at intermediate exciton-phonon couplings 
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unveils a striking non-monotonic behavior of D as a function of temperature, provided 
the lattice is actively modulating the hopping integrals in the TB Hamiltonian. This 
effect, reported here for the first time in the presence of an explicit environment, 
agrees with a long-known prediction made on the basis of a generalized master equation 
for the one-particle density matrix containing both dephasing and amplitude-affecting 
operators [43]. Intriguingly, we find that diffusive transport is faster at low temperatures 
with pure-diagonal noise [i.e. only on the site energies), but adding non-diagonal noise 
makes spreading faster at high temperatures. 

Importantly, our results on the diffusive regime at intermediate coupling flag a 
nontrivial interconnection between exciton spreading and collective (hydrodynamic) 
vibrational modes of the underlying lattice. More precisely, our findings strongly suggest 
that the observed non-monotonic behavior of the diffusion constant versus temperature 
is related to the presence of long-wavelength acoustic modes. This conclusion is 
reinforced by a comparison with the results of Troisi and Orlandi obtained in a similar 
semi-classical model with purely off-diagonal dynamical disorder [45] . In fact, they found 
a monotonic decrease of D with temperature in a model that lacks collective vibrational 
modes by construction, as their TB Hamiltonian is modulated by the dynamics of an 
ensemble of independent, disconnected classical oscillators. In fact, the signatures of 
long-wavelength hydrodynamic modes are clearly recognizable in the equilibrium power 
spectra S(k, to) of the exciton-coupled lattice, as shown in Figs. 1. This strongly suggests 
that coupling to an extended dynamical system, such as the one we consider here, 
might be a necessary condition to obtain non-monotonic transport with temperature, 
in agreement with [43]. 

Our results on the role of the lattice in the spreading properties of a quantum 
excitation show that the presence of non-diagonal coupling is sufficient to suppress 
Zeno-like localization at high temperature. To shed further light into this phenomenon, 
we then computed a measure of quantum efficiency for different choices of the chain 
parameters. Our results clearly confirm that, when the hopping rates in the TB 
Hamiltonian are explicitly modulated by the lattice dynamics, the transport efficiency 
is no longer quenched at high temperature, as observed by some authors in the absence 
of an explicit environment [19, 20]. Moreover, we fold that nonlinearity in the lattice 
dynamics exerts a powerful boosting action on the efficiency at high temperatures, 
confirming recent results on the importance of spatial and dynamical correlation patterns 
within the noise bath [83]. 

Overall, the results presented in this paper allow us to conclude that the properties 
of excitouic systems interacting with noisy environments are subtly shaped by the 
specific properties of the exciton-phonon coupling and by the nature of the spatio- 
temporal dynamical correlations that characterize the underlying lattice. It would be 
extremely interesting to extend the formalism presented here to more complex systems, 
such as light-harvesting complexes, which are widely studied in the context of quantum 
biology [84, 85]. 
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